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Abstract 


Differential cross sections of deep inelastic scattering of charged leptons from hadrons are 
investigated by using the gauge/string duality. We consider vector mesons derived from 
different holographic dual models obtaining a general expression. We focus on the strongly 
coupled regime of dual gauge theories for different values of the Bjorken parameter. We find 
new predictions which are particularly interesting for differential scattering cross sections 
of polarized leptons scattered off polarized vector mesons. We also carry out a detailed 
comparison of the moments of the structure functions with lattice QCD results. 
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1 Introduction 


In a recent series of papers [1, 2, 3] we have studied deep inelastic scattering (DIS) of charged 
leptons from scalar mesons and from polarized vector mesons at strong coupling by consider¬ 
ing different holographic dual models based on flavor Dp-branes in the probe approximation. 
In those papers we have studied the planar limit of dual gauge theories corresponding to 
D3D7-, D4D6D6- and D4D8D8-brane models ([4], [5] and [6], respectively) 4 . In [1, 2, 3] we 
have investigated the hadronic tensor and, for polarized vector mesons, we have obtained 
the corresponding eight structure functions. We have derived general relations, i.e. model- 
independent relations, of the Callan-Gross type between different structure functions. This 
is very interesting because from that one can infer a sort of universal behavior which would 
also be expected for QCD in the planar limit. This universal property holds because, for all 
the probe flavor Dp-brane models we have considered, mesons are described as fluctuations 
of flavor Dp-branes in terms of the corresponding Dirac-Born-Infold action. In particular, in 

[1] we have studied DIS from dynamical holographic scalar and polarized vector mesons with 
one flavor, by considering the flavored Dp-brane models of references [4, 5, 6]. Moreover, in 

[2] we have carried out a non-trivial extension of these results for multi-flavored mesons in 
the planar limit and with the condition Nf <C N, where N denotes the number of colors and 
Nf the number of flavors. In that paper we also have worked out the corresponding next-to- 
leading order Lagrangians in the 1/N and Nf/N expansions. In [1, 2] we have focused on the 
Bjorken parameter regime where l/\/A <C x < 1, and therefore the calculations have been 
done by using the gauge/string duality in terms of the type IIA and type IIB supergravity 
dual descriptions of the gauge theories in the large N limit and at strong’t Hooft coupling 
A. 

In a more recent work [3] we have studied DIS of charged leptons from hadrons in the 
C i < 1/VA and x ~ e -v ^ regimes. In [3] we have focused on single-flavored scalar 
and vector mesons in the planar limit. This has been investigated in terms of different 
holographic dual models with flavor Dp-branes in type IIA and type IIB superstring theories 
[4, 5, 6]. We have calculated the hadronic tensor and the structure functions for scalar and 
polarized vector mesons. In particular, for polarized vector mesons we have obtained the 
eight structure functions at small values of the Bjorken parameter. 

In the present work we further extend our previous results by studying DIS cross sections 
of charged leptons from polarized vector mesons and carrying out a qualitative comparison 
between moments of structure functions calculated by using gauge/string duality methods 
and available lattice QCD data. We emphasize that the character of the comparison which 
we perform is qualitative since these structure functions have been obtained in the context 
of the large -N and large-A limits of confining gauge theories derived from their string theory 

4 The calculations carried out in [1, 2, 3] are inspired in the methods developed in reference [7], where 
DIS of charged leptons from glueballs in the planar limit of N = 1* SYM theory has been considered [8]. 
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dual models, which in fact do not lead to real QCD. At this point we can briefly describe some 
important differences between QCD and the holographic dual models we consider. We can 
start from the D3D7-brane model, which is the holographic dual description corresponding 
to the planar limit of the strongly coupled SU(N) AT = 2 supersymmetric Yang-Mills theory 
with fundamental quarks [4], Obviously, this is not the large N limit of QCD but it is 
related to it. On the other hand, we study the D4D8D8-brane model developed in [6] which 
at low energy is expected to be within the same universality class as the planar limit of 
QCD. At high energy this model behaves in a different way compared with QCD, since in 
that regime the Kaluza-Klein modes generated by the compatification of the x 4 -direction 
of the N D4-branes along the S 1 become relevant. Therefore these additional degrees of 
freedom, which are obviously absent in QCD, play a role in the dynamics of this dual 
gauge theory at high energy. In addition, this model [6] enjoys an SO( 5) global symmetry, 
which is not present in QCD. Thus, one can say that these Dp-brane models do not lead 
to holographic dual gauge theories with exactly the same properties as QCD in the planar 
limit. However, these holographic dual models represent infrared confining gauge theories 
and share dynamical properties with QCD at low energy. For instance, the Sakai-Sugimoto 
model realizes spontaneous chiral symmetry breaking. Also, the chiral Lagrangian has been 
derived from this model [6]. 

We also ought to mention that the comparison we carry out between our results and 
phenomenology unfortunately has the disadvantage that our calculations and the existing 
phenomenology and lattice QCD calculations belong to different parametric regimes. On 
the one hand, in our previous papers [1, 2, 3] we have done calculations within the Regge 
region. The asymptotic Regge region is relevant since it dominates total cross sections and 
differential cross sections when considering small angle. However, in the asymptotic Regge 
regime it is not possible to perform quantum field theory (QFT) calculations, neither one can 
use lattice QCD methods in order to obtain QCD amplitudes [9]. From this discussion one 
infers that it is difficult to make any quantitative comparison of these models and SU( 3) QCD 
data. In this way our aim is to compare qualitatively certain aspects of the phenomenology 
derived from these dual models with QCD. At the end, by taking into account all these 
caveats, we would like to see how our results compare with lattice QCD data. Indeed, this 
kind of comparison has proved to be qualitatively (and even in certain cases quantitatively) 
interesting in several situations as we shall explain in the last section of this paper. In fact, 
from a theoretical perspective it has given some important insights about the behavior of 
gauge theories in general, in the strong coupling regime where perturbative QFT methods 
cannot be used. This assertion is indeed very important because given a certain solution of 
supergravity and/or superstring theory, in principle, the gauge/string duality can be used 
to define a strongly coupled holographic dual SU(N) gauge theory in the large N limit. 
Suppose in addition that such a QFT has composite states. A natural way to study these 
states is in terms of the DIS calculations as described in [7, 9, 1, 2, 3]. Therefore, the interest 
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of these methods goes beyond the large N limit of QCD itself, i.e. it is useful for any other 
confining gauge theory with a string theory holographic dual model. 

This paper is organized as follows. In Section 2 we introduce and derive a number of 
relevant expressions necessary to construct the DIS cross sections of charged leptons from 
hadrons of spin s = 0, | and 1 in a generic form. This is a fully general expression for the DIS 
differential cross section which holds for any infrared confining gauge theory at any value 
of the rank of the gauge group N (i.e. beyond the planar limit), and at any value of the 
coupling. In Section 3 we focus on the case of polarized vector mesons and we present these 
derivations in detail. Thus, in this section we derive the DIS cross section from leptons off 
vector mesons which are unpolarized, longitudinally polarized, transversally polarized and 
with partial polarization. We also study the corresponding helieity amplitudes of the forward 
Compton scattering directly related to the DIS process and discuss on physical implications. 

In Section 4 we derive general (i.e. model independent) expressions for the DIS cross 
sections at large ’t Hooft coupling, in the three mentioned regimes of the Bjorken parameter. 
In general terms the hadronic tensor can be derived from the operator product expansion 
(OPE) of two electromagnetic currents inside the hadron. At weak ’t Hooft coupling the 
OPE is dominated by approximate twist-two operators, and it corresponds to the scattering 
of a lepton from a weakly interacting parton. However, in this section we consider the 
strong coupling limit. In this case the OPE which leads to the hadronic tensor is dominated 
by double-trace operators, and it corresponds to the scattering of a lepton from the entire 
hadron. Within this strongly coupled regime we consider the parametric range l/\/A 

1 where supergravity describes DIS. Furthermore, we consider the case where e _v/ * C i C 
i/Vx, in which excited strings are produced. Also, we study the case for x ~ e , where 
the produced excited strings have a size comparable to the scale of the AdS space. 

In Section 5 we discuss our results derived in the previous sections such as differential 
cross sections in terms of the Bjorken parameter and the fractional energy loss of the lepton, 
as well as the structure functions and their moments, and carry out a comparison with 
moments of structure functions for the pion and the p meson obtained from lattice QCD. It 
turns out that our results of the first three moments of the F 2 -structure function of the pion 
and the iq-structure function of the p meson are in very good quantitative agreement with 
the corresponding ones from lattice QCD. 

2 Deep inelastic scattering cross section 

In this section we first introduce some definitions of kinematical variables which will be useful 
for the rest of the work. Then, we introduce the leptonic and hadronic tensors and give an 
expression for the DIS differential cross section which holds for any value of the coupling. 

The DIS process is schematically represented in figure 1. 
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Figure 1: Schematic representation of DIS. A lepton with momentum k interacts with a hadron 
with momentum P through the exchange of a virtual photon with momentum q. 

We consider the following definitions which are the same as in reference [10] 5 

• M is the mass of the hadron, being the on-shell condition M 2 = — P 2 . 

• E is the energy of the incident lepton. 

• k is the momentum of the incident lepton k 11 = (E, 0, 0, E). 

• E' is the energy of the scattered lepton. 

• k' is the momentum of the scattered lepton k = E'(l,sm6 cos </>, sin 9 sin0, cos 9). 

• P is the momentum of the hadron = (M, 0, 0, 0) in the rest frame 6 . 

• q is the momentum transfer or momentum of the virtual photon q = k — k'. 

• v is the energy loss of the lepton v = E — E' = —-jy. 

• y is the fractional energy loss of the lepton y — — — 

2 

• x is Bjorken parameter x = 

• t is a parameter defined as t — 

In what follows we have assumed that the leptons are almost massless. 

The DIS differential cross section can be written in terms of the so-called leptonic and 
hadronic tensors and W respectively. Thus, we have 


dx dy dcf> 



( 1 ) 


5 Notice however that we use the mostly-plus metric =diag(—, +, +, +). 

6 Fixed-target experiment. 
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where e denotes the electron charge. The leptonic tensor depends on the incident lepton 
beam and it can be easily derived from perturbative QED, leading to 


C = [2(k»k' v + k v k*)-2rTk-k , ] + [-2i(? , ' af) q a sif i ] , m z 2 ~ 0 (2) 

2 

_ 1W i 

l sym ' L ant ? 

where and denote the symmetric and antisymmetric parts, respectively, sf is the 
lepton spin four-vector defined by boosting (0, s{) from the rest frame, or equivalently by 
using the spinor states u and u which leads to 

Si = u(k, si)Yl5u(k, si) • (3) 

On the other hand, the hadronic tensor describes the internal structure of the hadron. It 
cannot be calculated by using perturbative QCD because there are soft processes involved. In 
this paper we focus on spin-one hadrons since they have a very rich structure given in terms 
of eight structure functions which follow from the most general Lorentz-tensor decomposition 
of the hadronic tensor. We also write the results of the structure functions for the case of 
scalar mesons since we shall discuss them in section 5. 

In general terms, the hadronic tensors for a spin s — 1, \ and 0 can be written as follows 


K = w%r + w™ + w£, (4) 

K-i = wp+w*, (5) 

K = K F - ( 6 ) 

where 

K F = - ~^P.P V , (7) 

Wpv = biTfj,, v + U^ v ) — U^y) ~^( s nv ~ tfj,y) , ( 8 ) 

m 'is = -jp-^x^s°-j^^^q x (P-qs°-s-qP°)- ( 9 ) 

The superscripts refer to the symmetry properties under the exchange of Lorentz indices (S 
and A) and to the usual name given to the structure functions involved in each term (TVs, 
Vs and Vs). In addition, we have used the following definitions [11] 

>> = -Cq-C-^(P- q) 2 k)Vhv , ( 10 ) 

s pp = , p 2 ■ C q ■ ( - \{p -(i) 2 k)PhPvi (h) 
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V = WpTtf - |( p ' «ft p - + ((« ' OftC. + Pft) + c.c.)' 

(12) 

U ^~ P . q [(CC + CU 3 3 PpPv] , 

(13) 

s ° = up, . 

(14) 

K = 1 — 4 x 2 t , 

(15) 

where is the polarization vector of the hadron. Then, one can write a general expression 

for the differential cross section of the DIS of a charged lepton from a hadron in terms of 
the symmetric and antisymmetric parts of these tensors 

72 4 

a — e v \//^ w SF i /^ w Aa i y w w Sb ] 

dxdyd(f) 16ttV y [ sym ^ ant ^ sym 

(16) 


where <f> is the azimuthal angle. The differential cross section depends trivially on </> but 
non-trivially on x and y, which can be chosen as two independent parameters describing the 
DIS process, together with the spin polarizations of the lepton and the hadron. It is worth 
to note that in the case of spin-zero hadrons only the first term in the square bracket of 
Eq.(16) contributes, while for the spin-1/2 case the second term also appears. The last term 
contributes for spin-one hadrons and contains four additional 6* structure functions. 


3 DIS differential cross sections and structure func¬ 
tions 

In this section we derive the DIS differential cross sections of charged leptons scattered from 
hadrons with spin s = 0, ^ and 1. We consider first the general case, and then we study 
different polarization states of the hadron beam. 


3.1 The general case 

Let us begin with the first term of Eq.(16), which is present in hadrons with spin 0, | and 
1. We obtain the following term expressed in terms of the F l and F 2 structure functions 


Vi y xy 2 F 1 +[l-y + x 2 y 2 t\ F 2 


167 r 2 q 


47r 2 g 4 
e 4 M E 


| xy 2 F ] + [1 — y] F 2 |, 
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where in the last step we have approximated |t| ~ 0.' The second term, which appears for 
both s = \ and s — 1 spins is given by 


167r 2 q 4 


»C,w* = 


e 4 ME r 2 

91 


2 x ( Sh ‘ s i) + (q ■ s h)(q ■ si) 


+2xy 2 g 2 


(p-q) (p-q) 2 

(s h -si) ( q-s h )(p-s t ) 


L (p-q) 


(p ■ qf 


(18) 


This term depends on the structure functions g\ and g 2 . Notice that in the expression for 
the first term we contract only symmetric factors, while in the expression for the second one 
only antisymmetric tensors are contracted. In addition, in the last calculation we have used 
the identity e„ uXa = -2 (rj^rjP - C'^a)- 

Let us now calculate the most general form of the third term which only appears for spin- 
one hadrons, and can be written in terms of the 6* structure functions. Firstly, we need to 
calculate the contractions 


//W' r lV” Q + . lfJ-v 

l sym 1 5 l sym °5 (j sym L Hv 5 L sym Li 


(19) 


In order to simplify the notation we define a pseudo-scalar product (••) given by 

/, , \ _ 3 (h ■ (*)(&2 • C) + c - c - 

{hk2) = 2 - (p^r -’ (20) 

where c.c. indicates complex conjugate. This pseudo-scalar product is linear on both sides 
and it is also commutative. In addition, since we take the \t\ <C 1 limit we set k — 1. Then, 
we obtain 

1 / .1 

3' 

4 

’3* 

4 


r 

' J sym‘ yv 


[2 («" + k"k'“) - 2 rTk ■ k')- jEm(i q?) 

;(k ■ k') l(qq) - 1] 


-MExy[(qq) - 1] , 


( 21 ) 


p iU o 
l sym s 


(2 (kV + k-k'“) - 2 rTk ■ k '\-prr^s('J ' C 9 ■ C - • q?)p^. 


4 \2(P ■ k)(P ■ k') _ P 2 (k ■ k') 

3 [ (P-q) (P-q) 

1 


i(qq) - 1] 


--ME [(qq) - 1] 


L v 


— 1 + xyt 


( 22 ) 


t 


= [2 (¥ L k' v + k v k ,fl ) - 2 rf'k ■ k r 


2 (P-q) 2 


x 


'This is the equation (5.2) of reference [10]. An extra factor 27 t appears in the denominator because we 
are calculating as in [10]. In reference [11], on the other hand, it has been calculated = 27T d ^ d y d( p ■ 



' - \{P • q)P»Pv + ((<? ' C)(^C + P„Q + c.c.y 

— / \ _ ^sym S ■ P ) + ~(qk)(k' ■ P ) 

= 2 Q - 1 + xV^ - (g£/) - (gfc)(l - y) , 

id 

= \ 2 (k“k"'+k-n - *r* • fc'i-2— [<c;c+co - §«v -Iw 


W. g) - 2 M 


P" r _ 

^ C'll'KKJ ' U.U 


-P"s„ 


(P-g) \ (P ■ q)[(qq) ~ l} sym ^ (qq) - 1 sym ^ 


= -ME —y(kk') + (-1 ) + 2 x 2 yt 


Finally, by writing the cross section contribution from the 6 * structure functions as a\b[ + 
« 2&2 + « 3&3 + CI 4&4 one hnds the following expressions for the ag coefficients 

e 4 MEe 4 

ai ~ I 67 r 2 g 2 ^ l l ‘y"> r i"' ~ 2 2 3Ti 2 q 4 ,r -V — ’ (^5) 

6 ^ y 

a 2 = — 167 T 2 g 4 6 ^ ym ^ S>1U ^ 

MEe 4 
2 3 3 2 7r 2 q 4 

[2 ((qq) - 3) (1 - j/) + j/(gfc') + y(l - y)(gfc) + 2y 2 (kk') - 2 (4 - (gg» xVt] , 

(26) 

a , =--- y. p" ( s -u ) 

3 167r 2 g 4 2 sym V yu yu 


MEe 4 
2 2 3ttV 


(1 - j/) (gg> - t/ 2 (A:A;') + xVi ((gg) + 1 )] , 


6^ y 

a4 = ~ 167 r 2 g 4 2 ~~ ^ 


MEe 4 
2 3 37r 2 g 4 


[2 (1 - y) (qq) ~ y(qk') - (gfc)j/(l - y) + 2 (qq)x 2 y 2 t 


By using the linearity of (••) we obtain 


- V P" IF 56 = 

167r 2 g 4 y s v m w 

MEe 4 f, 0r , , 


127r 2 g 4 


6 iX|/ 2 [(gg) — 1 ] + — [(2 — 3y)(gg) + y(2 - 3 y)(qk) + 2 g 2 (/c/c) - 6(1 - y)\ b 2 
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+ [(1 — y) (qq) ~y 2 (kk) +y 2 (qk)} b 3 + i [(2 — jr) (qq) +y(-2 + y)(qk)]b i + 0(t)X . 

(29) 

Now, in the yf—t <C 1 limit we can relate the values of (qq), ( qk ), and (kk), since in this 
limit '-jjo — where the Lorentz index p runs from 0 to 3. Then, in the hadron rest frame 
we obtain the relations 

q^q u \(q^ l k v + k fl q u ) E L k u 

M 2 (q 0 ) 2 ~ M 2 q°E “ M 2 E 2 ’ 


2 


qtiqu 

(p 7 # 


C// C/y 


(<fk v + E L q u ) _ 0 k^k u 
(P ■ q)(P ■ ky^ v (P-k ) 2 ^ l/| 


3 (g-C)(g-0 + c.c. = (P • g) 3 (g • C)(fc • C) + C.c. (P • g ) 2 3 (k • C)(fc • 0 + c.c. 
2 (P-g ) 2 (P ■ k) 2 (P-g ) 2 (P-k) 2 2 (P-q) 2 

and 


(gg) = j/(gfc) = ?/ 2 (*P) • 


Thus, for <C 1 the part of the cross section associated with the 6 j structure functions 
of the hadronic tensor reduces to a simpler form 

e 4 M F,p a 

16^4 V - 1) [W + *2(1 - »)] + ® (v^) . (30) 

Notice that in this limit the contributions to the cross section coming from 63 and 64 are 
sub-leading in t even if & 3 and 64 are not necessarily sub-leading themselves. Thus, the full 
DIS differential cross section from spin-one hadrons to this order becomes 


da 


spinl 


dx dy d(j) 


| [ x v 2f i + (! - y) F 2] + ^((gg) -1) [xy 2 bi + (1 - y)b 2 ] 


+y 2 gi 


MEe 4 
47 r 2 g 4 


+y 2 gi 


2 x ^—LL + (V ' s h)((l ■ s i) 


(p-q) 


(p ■ qf 


xy 


Pi + x((w) - X ) & i 


2x ( Sh ' Sl ) | ( q ' s h)(y- s i) 


(p-q) 


(p ■ qf 


+ 2 xy 2 g 2 

+ (1 -y) 

+ 2 xy 2 g 2 


(sh-si) (q-s h )(p-si) 


(p-q) 


(p ■ qf 


f 2 + -((gg) - 1)62 

(q-s h )(p-Si) 


L (p-g) 


(p ■ qf 


(31) 
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3.2 Hadron polarizations 

Now, we calculate the dimensionless factor 


3 (g-C*)(g‘0 + c-c- = o c r 

2 (P • q) 2 (P • q) 2 ^ ’ 

for different polarizations of the hadron in order to obtain the DIS differential cross section 
(h . c ^ d( p for the particular cases which may be phenomenologically relevant, namely: when 
the hadron beam is unpolarized, longitudinally polarized, transversally polarized and also 
we study partial polarizations of the hadron beam. In all cases the longitudinal or 5-axis is 
defined by the spatial components of k, in other words, the direction of the incident lepton 
beam. 




Unpolarized hadron beam 

In order to describe unpolarized hadron beams we have to average over polarizations obtain¬ 
ing 

TO = + P»Pu) ■ (33) 

Then, 

(< qq)un P oi = Cm) = + (P ■ q) 2 } = i - 4 x 2 t = K ~ 1, (34) 

thus, we obtain a null contribution from the bj structure functions as expected for an un¬ 
polarized target [11]. In particular, for a longitudinally polarized leptonic beam we have 
Si = Hik with Hi = ± 1 for positive and negative helicity, respectively. Then, we obtain 

da MEe 4 r 9 

dxdydi ) = 4^V ^ 1+ 2 ^‘ 


Longitudinal polarization 


Now, we consider the case when the hadron beam is polarized in the direction of the incident 
beam and call it 5-axis. Thus, we can take the polarization to be 


Cl 


M 

7! 


( 0 , 1 ,^, 0 ) 


->■ s a 


H h Mz , 


(35) 


where Hh = ±1 indicates that the hadron polarization is parallel or anti-parallel to the beam 
direction, respectively. Then, we obtain the following expressions 


(<? ' Cl) 


ME' sin# 
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(q-Cl)(q-Cl) 

(qq)LP 


M 2 E ,2 sm 2 9 q 2 M 2 


3 


2 

1 


2 (P-q) 2 


2 

M 2 E ,2 sm 2 9 


[(1 -y) + tx 2 y 2 ], 

= — 6x 2 t[(l — y) + tx 2 y 2 } 




0, 


(36) 


In this case the factor ( qq ) — 1 only gives a minus sign. In particular, for a longitudinally 
polarized lepton beam we obtain 


da 

dx dy d(j) 


MEe 4 
47r 2 g 4 



+ (!-?/) 



HiH h {2 - y)yxg\ + C?(i)j . 

(37) 


Transversal polarization 


Since the system under consideration has azimuthal symmetry we can choose any polarization 
for the hadron in the (x, y)-plane. For transversal polarization we choose for example the 
x-axis and set 


M 

Ct — iH h ) —» s a = H h Mx 


H h = ±1. 


Then, we obtain 


(q • C t) — 


(1 -y) 


-t 


M 2 

7f 


M 2 r l(l-y) . , iH h 

Sill 0 


— x 2 y 2 sin0 + iH h (-f- xy 


2 xt 


(9-Ct)(9-Ct) = 
(qq)rp = 


V2[\ ~t 

M 4 r 1 
2t 2 _4x 2 
3 1 M 4 r 1 


2xt 


2 (P-q) 2 t 2 


— t (1 — y)sin 2 (/> 

J 

— t(l — y)sin 2 (/> 


4x 2 


6a;2 1 4^2 - i ( 1 _ y) s[n2< i 


3 

2 


(38) 


(39) 


and the dimensionless factor becomes (qq) — 1 « f. In particular, for a longitudinally 
polarized lepton beam we obtain 


da 


MEe 4 


dx dy d(j) 


xy 


F\ + -&i 

o 


+ (i-y) 


F2 + ~h 

o 


47r 2 g 4 

PHiH h 2x 2 \fM\J\ - y(yg x + 2y 2 ) cos^j . 


(40) 
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Partial polarizations 

The b t structure function terms in Eq.(30) and in particular the dimensionless factor (qq) — 1 
can be written in a more general form 

ypr rp 4 _ 

W v l nm W H£ = Tr(p ■ As) [W + Ml -</)]+ O (V=t) . (41) 

where As = diag(l, 1, —2) and p is the spin-one density matrix that accounts for the 
possibility of having a statistical production of the hadron beam. Notice that 

Tr(p - As) = [(Si) + (5 2 ) - 2(Sl) + 3(5,)] . (42) 

We can study particular cases with partial longitudinal and transversal polarizations. In the 
A-axis basis we have p U n P oi = (1/3) diag(l, 1,1) and 


PpLP — ^ 1 2(1 — Hp) 0 ; PpTP — 


2 - Hi 0 -2 + 3 Hi 


-2 + 3 Hi 0 2 - Hi 


where Hi and Hi are the fraction of longitudinally and transversally polarized hadrons 8 , 
respectively. Here, H 2 = 1 represents the totally polarized case and H 2 = | describes an 
unpolarized hadron beam. In the former notation we identify 

{qq) P LP = -3 0 -Hl'j (qq) L p + 3(1 -Hi) (qq) unpol , (44) 

( qq) P TP = -3 Q - H^j ( qq) TP + 3 (l - H%) ( qq) U n P oi , (45) 

where the subindex p stands for partial polarization. Then, the results for this part of the 
DIS differential cross sections corresponding to a partially polarized target are 


MEe 4 / 2 


167r 2 g 


- v P" W bb = - 

4 y u sym v ' pis 2 2 7l 2 q 4 


[W+Ml-sO] +0(v / =i). (46) 


for partial longitudinal polarization, and 


_T_ y r w » = 

16irV 9 2%V 

for partial transversal polarization. 


^Hi — [b±xy 2 + b 2 ( 1 — y)\ + O (V—t) , (47) 


8 Again, we are writing transversal polarization along the £-axis. However, by choosing the y-axis instead 
it would lead to the same results. 
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These expressions lead to the final results of DIS differential cross section for partially 
polarized hadrons as a generalization of the previous cases. For a partially longitudinally 
polarized target Hl we obtain 


da 


MEe 4 


dx dy dcj) 


xy 


F '+ I 3 


47r 2 g 4 
-H t H h (2 - y)yxg 1 }, 


Hi h 


+ 0--V) 


--Ilf \b 2 


(48) 


while for a partially transversally polarized target Ht we obtain 


da 


dx dy d(f) 


MEe 4 
47r 2 g 4 


xy 


F\ 


-.--Hi) h 


+ (i -y) 


+HiH h 2x 2 \fM^/\ - y(yg x + 2g 2 ) cos qb 


1 

2 


- - -HI I b 2 


(49) 


We have assumed longitudinally polarized leptonic beams Hi in both cases. Once again, by 
assuming that F 2 (x) = 2xF { (x) and b 2 (x) = 2xb\(x) are satisfied, we recover the formulas 
shown in [11], 


3.3 Helicity amplitudes 

Deep inelastic scattering is related to the forward Compton scattering. Let us briefly describe 
the forward Compton scattering by using the virtual photon momentum q and the hadron 
momentum p in analogy with our DIS notation, and the variables A and X' which denote 
the initial and final hadronic hclicities. In the forward Compton scattering both the initial 
and final states of the system contain a photon and a hadron. The related amplitude is 
a function of a tensor T /i!y defined in terms of the matrix elements of two electromagnetic 
currents by the formula 

(T„„) u , = i J (fW’-'te A'l T (J„(x)J„( 0)) |p. A), (50) 

which has a structure similar to the hadronic tensor W^ u involved in DIS. On the other 
hand, the definition of W^ v in terms of the currents is given by 

(WV) M , = 2 - 1 A'l I J„(x), J„(0)] |p. A). (51) 

Since both T^ u and W /IU have the same Lorentz symmetry transformation properties, they 
have the same tensor structure decomposition. Thus, the information they contain can be 
encoded in their structure functions in the same way. 

DIS and forward Compton scattering are related to each other by the optical theorem 
in quantum field theory as it is schematically shown in figure 2. It implies that twice the 
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imaginary part of the forward Compton scattering amplitude leads to the total DIS ampli¬ 
tude. Therefore, the imaginary part of the structure functions give the W /w structure 
functions times an extra numerical factor. Actually, this relation between the two processes 
is one of the basis of our construction to obtain the structure functions. As we will see in the 
next section, it is possible to describe a holographic dual version of the forward Compton 
scattering from the point of view of supergravity and string theory via the gauge/string 
duality. 



Figure 2: The relation between DIS and forward Compton scattering processes as given by the 
optical theorem in terms of Feynman diagrams. 


The Compton amplitude can be analyzed for different helicities of the incoming and out¬ 
going particles. The amplitudes in each case are commonly denoted by A hH . h i H i, where 
the subindex h corresponds to the helicity of the photon and H labels the helicity of the 
hadron. Primed variables indicate helicities of final states, otherwise labels indicate helicities 
of initial states. These helicities can take three different values in the spin-one case, —1, +1 
and 0, and they must satisfy the condition h + H — b! + H'. As shown in [10, 11] they can 
be calculated as 

Ah.H,h'H 1 = 4* e h'Tin>{s h) , (52) 

where the e’s are the photon polarization vectors and Sh is the hadron spin. For spin-one 
hadrons there are eight independent helicity amplitudes. By considering a general kinemat- 
ical regime, and by using the notation from [11], where <23 = b- 2 /3 — b 3 and <24 = b 2 /3 — 64 , 
these helicity amplitudes can be written in the following form 


A 


++,++ 


A 


+o,+o 


A 


+0,0+ 


^ (1 — 4:X 2 t)bi xt A 0 

F\ - y - -— - — a 3 - 9l - Ax 2 tg 2 

^ 2(1 — Ax 2 t)bi 2 xt 

F - + 3 


~ F 1 


2 64 

IT’ 


^[2*( gl +ri + l+ (i - 4 fb 


(53) 


(54) 
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(55) 


- 4 +-, 00 


^4+-,-+ 


A 


0+,0+ 


A 


00,00 


[2x(^i + 32 ) + y + y ] , 


^+-,+- = *1 


(1 — 4x 2 t)6i xt a n 

---y as + 9i + 4x 2 ^ 2 


F &i , 

- ^1 - y + 9i, 

= V-< [2x(0i + 02 ) - y-|-] 

- + 92 ) - y - y] , 

= —2xta 3 ~ 0 , 

(1 — 4a: 2 f)F 2 (1 — 4x 2 t)6i (3 — 12o: 2 t + 16x A t 2 )b 2 

- —r\ H---h 


2x 

A(—2x 2 t)bs 2xt. 2 

+ 3x +y-(l~ 4ar 2 t)6 4 


18a; 


~ — Fi 


£2 , _ £l 

2x 3 6a: 


(1 — 4a: 2 t)F, 2(1—4a: 2 f)&i (3 — 12a: 2 t + 16a: 4 t 2 )5 2 

_Fl +-bT-3- + -(w- 

JAzAAh _ « ( i _ 4 *^ 

3a: 3 V ; 

F 2 2b\ b 2 

~ Fl+ 2i~~ + 3V 


(56) 


(57) 

(58) 


(59) 


(60) 


Notice that the first line in each equation (53)-(60) coincides with the corresponding one 
in equations (7) in reference [ll] 9 , and we write them here for completeness. On the other 
hand, the second line in each of these equations corresponds to the limit |t| <C 1. 

These amplitudes are useful tools in order to study the hadron structure, specially in the 
cases with s > 1/2. They provide a simple way to understand the meaning of the structure 
functions in theoretical terms, and give important insights about how they can affect the 
scattering amplitudes and how to measure them. This is the reason why we shall explicitly 
show the form of these hclicity amplitudes in terms of the structure functions we found by 
using the gauge/string duality. 

For instance, in the spin-1/2 case one has A + ^ j+ ^ ~ F\ — <?i and A + ± t+ ± ~ F\ + g 1 , which 
indicates that iq is the cross section for a transverse photon scattered off an unpolarized 
target, while g\ is the spin asymmetry in the scattering cross section for a transverse photon 
[10]. One can also say that g\ +g 2 is proportional to the single hclicity flip amplitude ~ A+^ot- 
For spin-one hadrons a similar analysis can be carried out. The function F\ ^4++,++ + 
7L +0i+ o + A + _ )+ _ can be interpreted in the same way. Also b\ ~ A ++i++ + A + _ i+ _ — 2A + o !+0 

9 This is easy to see by writing k in terms of x and t, by taking into account an overall sign in the definition 
of the four-dimensional metric, and by noting that the definition of v = P • q in reference [11] differs from 
the one used in the present work which is v = — P ■ q/M. 
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has an interpretation from a transverse photon scattered off a polarized spin-one target [11]. 

The combination gi+g 2 also appears in the A. + _ ; oo amplitude, and at first order in t this is the 

only place where r/ 2 is present. The combinations F\ — £>i/3 and -F 2 — 5 2 /3 define A 0 +,o+ and 

will be important when we present our structure functions. Also, the A + __ f ~ —& 3 + 6 2 /3 

is an interesting amplitude (even being sub-leading in t) since it characterizes the double 
hclicity flip. As pointed out in [10] 10 this quantity does not receive contributions from quark 
interactions at first order in a s because quarks cannot change their hclicity by two 11 . 

Finally, all these linear combinations of structure functions that define the hclicity ampli¬ 
tudes are positive, and this is a consistency check for the structure functions we obtained in 
our previous papers [1, 2, 3] in terms of the holographic dual models. 


4 DIS differential cross sections from the gauge/string 
duality 

The present approach is based on a method developed by Polchinski and Strassler in [7], 
where the structure functions corresponding to glueball states were derived in terms of the 
gauge/string duality. They consider the planar limit of the Af = 1* SYM theory obtained 
from a deformation of the Af = 4 SYM theory. The Af = 1* SYM theory has a running 
coupling and becomes conformal when the energies involved are higher than a color con¬ 
finement scale A. Therefore, this SYM theory and QCD share some important qualitative 
properties: they exhibit confinement of the color degrees of freedom and a mass gap in the 
infrared. In the context of the gauge/string duality it was shown that the Af = 1* SYM 
theory is dual to strings propagating in a certain deformed AdS§ x S 5 background [8]. In 
order to describe DIS, it has been calculated the scattering amplitude of a process where 
a particular graviton perturbation coming from the boundary couples to the dilaton, which 
is dual to a glueball state [7]. In addition, the optical theorem is used to relate DIS and 
forward Compton scattering. Therefore, in the gravity dual description the graviton and the 
dilaton perturbations are taken to describe both initial and final states leading to a process 
which turns out to be dual to the forward Compton scattering described in the previous 
section 12 . It is important to emphasize that in the large-x case the supergravity description 
is appropriate to derive the structure functions, while in the small -x regime the string theory 
description is unavoidable. 

The calculations of DIS from vector mesons have been presented in our papers [1, 2, 3]. 
We have obtained the leading order behavior of the eight structure functions involved in the 
hadronic tensor W^ u for a spin-one target, namely F\, F 2l bi, b 2 , & 3 , & 4 , g\ and g 2 in terms of 

10 In the mentioned paper this is the so-called A(:r) structure function. 

11 Recall that a s = 3qc£>/47t. 

12 For the large- 2 : case there is also a sum over intermediate dilaton states. 
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x , A) and some other parameters which define the holographic dual description. 

The Bjorken parameter a: is a very important kinematic variable, and its physical range 
0 < x < 1 can be divided in four regions 13 . In the region l/\/A <C x < 1 (which we call 
the large -x region or range A) the ten-dimensional Mandelstam variable associated with 
the intermediate state s ~ (a' 2 A) _1 / 2 (l/a; — 1) indicates that only massless string states 
are exchanged and thus we only have to consider the low energy limit of string theory, 
leading to a supergravity description. The second kinematic region (range B) is reached 
when e~ ^ <C x <C l/\/A. Thus, in this region the dynamics of string theory provides 
the appropriate description. Fortunately in this region, it is possible to think of locally 
interacting strings, in the sense that the characteristic length scales are smaller than the AdS 
radius, leading to some important simplifications [7]. This is the region we call small-x (but 
not exponentially small). The third region (range C) corresponding to x ~ e~ ^ is important 
because the non-locality of the string scattering is described by a diffusion operator. This 
operator leads to the convergence of the first moments of the structure functions [7, 9]. For 
smaller x values, there is fourth region which is not accessible with the present approach 
since it would require to consider strings propagating in a curved background. Finally, we 
point out that for actual experiments large x generally corresponds to 0.1 < x < 1, while 
the small x measurements are supposed to describe the 0 < x < 0.1 region. 

In the following sections we present our results for the structure functions in the first two 
regimes and show how they affect the differential cross section of the DIS process in the 
different polarization cases. We also consider the corresponding hclicity amplitudes. 


4.1 Range A: 1 /a/A «x<1 

The structure functions found in our previous paper [1] in this regime of the Bjorken param¬ 
eter at leading order in t are 

F\ = 

F 2 = 
h = 
b 2 = 

h = 
b A = 

13 Recall that x = 1 corresponds to elastic scattering. 


~ x ) » 


~ x ) > 


12a: 3 

1 




6a; 3 

1 


~ x ) > 


4a; 3 

1 


- 4x ) > 


2a; 3 

1 


24a; 3 

1 


(61) 


^ ( x ) 12^ ( ~ 1 + 4a;) ’ 
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Model / Parameter 

P 

a 

P 

D3D7 

7 

2 

-2 

D4D8D8 

8 

3/2 

-3/2 

D4D6D6 

6 

3/2 

-3/2 


Table 1: Parameters for each Dp-brane model. 


gi = aSF (x)-^(-7 + 6x) , 


t 


8x 2 

1 


where 


92 = A (a;) —(3-3a;), 

A SF (x ) = A% F n 2 p Q 2 a' 4 R 2p ~ 6 (^)\ T+2n+5 (l - a;) 7 " 1 , 


(62) 


is the only model-dependent factor and it is written in terms of q 2 and the different variables 
involved in the string theory calculation: the Dp-brane tension g p , the charge of the hadron 
Q, the confining scale A, the sphere radius R, the string constant ol and the exponents 

A 2 + £(£ + p- 5) A _ l-2a-f3(p-3)/2-p + 4 n a - j3 - 2 


7 2 = 


B = 


B 2 2 1 2 

related to the general induced asymptotic metric on the probe Dp-branes 


ds 2 = 


(^) rj^dx^dx" + [dr 2 + r 2 dtt 2 A 

Vii/ Vii/ 


(63) 


(64) 


Furthermore, we have defined n = Notice that l is the usual spherical harmonic index 
that appears in the solutions of the scalar and vector mesons, and A FF is a dimensionless 
normalization constant which depends on the model [1, 2], As explained in [2], for the three 
considered Dp-brane models the values of the parameters are shown in Table 1. Using the 
redefinition fiy = - A -^A SF (x)Wi we obtain 


Fi = (1 — x) 
b\ = 3(1 — x) 

h = ^(1 -4x) 
^ = ^( 3 - 3 ^) 


Fj = 2(1 — x ), 
b 2 = 6(1 - x) 

h = -(1 - 4x), 

~ 3 xt, „ 

9i = —(-7 + Qx) 


0. 


This means that in this regime we obtain the relations [1, 2] 


F 2 — 2 Fi , b 2 — 2bi , b\ — 3F X , g 2 — — iy , 6 4 — —2 b 3 

Ax 
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Note that the first two are typical Callan-Gross type relations, except for the lack of a factor 
x , which in this regime is understood as if the lepton is scattered by the entire hadron due 
to the strong coupling and the specific kinematic regime. By looking at the g 2 -dependence 
we obtain a power-law decay, which is not exactly the same as in perturbative QCD. This 
was already observed in [7] where it was pointed out that at weak ’t Hooft coupling the 
leading contribution comes from twist-2 operators, while at strong coupling the dominant 
contribution is given by double-trace operators. 

For completeness we also write the results for the scalar mesons, which as mentioned only 
have two structure functions 

/ A 2 \ 7+1 

F, = 0 , F -2 = A s 0 calar fi 2 p Q 2 a'*R 2p “ 6 l — \ o; 7+3+2n (l - x)\ (65) 

where A s 0 calar = 2 27+4 £? 27+2 7t 5 | q | 2 1 c x | 2 » while the constants c t and c x have been 

dehned in [1], 

By plugging the structure functions from the holographic vector mesons in the general 
spin-one DIS differential cross section we find 


da 


dx dy dxj) 


MEe 4 
47r 2 g 4 
MEe 4 


(qq) [xy 2 + 2(1-7/)] F 1 + 2 xy 2 g 2 


( s h -si ) (q-s h )(P-Si) 


Fi \ (qq) [xy 2 + 2(1 - y)] + -y 2 


47r 2 g 4 
MEe 4 A SF (x)( 1-x 


L (P-q) ( P-q ) 2 

( s h -si) (q ■ s h )(P ■ si)- 


L (P-q) 


(p ■ qy 


4n 2 q 4 


12a; 3 


x 


(qq) [xy 2 + 2(1 - y)] + -y 2 


(sh-si) (q-s h )(P-Si) 


(P-q) 


(P-qf 


( 66 ) 


For a longitudinally polarized lepton beam where s; = Efik with Hi = ±1 for positive and 
negative helicity, we can rewrite this expression in the following form 


da 

dx dy d(j) 


MEe 4 A SF ( x) 

487 r 2 q 4 x 3 X 


| (qq) [xy 2 + 2(1 - y)] 


m 

2 ME 


*h (V k n 



.(67) 


We can study the behavior of this result for each polarization as follows. 


Unpolarized target: (qq) unpo i = 1, s h = 0. 

da | MEe 4 A SF (x) 


dxdydcf) 


unpol 487 f 2 q 4 x 3 


(1 -x) [xy 2 + 2(1 - y) . 


( 68 ) 


• Longitudinally polarized target: ( qq)LP = 0 , Sh = (0 ,H h Mz) (and using cos 9 ~ 1). 


da 

dxdydcf) 


LP 


0. 


(69) 
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This is the result at leading order. It is easy to see why this happens: since ( qq )lp = 0, 
the bi = 3 F{ relations for i — 1,2 imply that the b-F part of the cross section vanishes. 
Now, as pointed out in [10] Sh is dotted with k or q both of which have similar 0 and 
3 components at leading order. Thus, we replace it by —H^p 14 without changing the 
results. This implies that the ^-piece becomes sub-leading in this case, and as g\ ~ 0 
at first order, therefore the cross section vanishes. 


• Transversally polarized target: ( qq)rp = 3/2, Sh = (0 ,H h Mx) 


da 

dx dy dcj) 


TP 


MEe*A SF (x) (3 

48t rVx 3 1 ; \2 


xy 2 + 2(1 - y ) 


0////4 M 
2 \q\ 


(1 — y) cos (ft 


(70) 


Note that the M/\q\ = \f—t factor makes the spin-dependent term sub-leading. The 
reader could think that it is not a reliable result since in early steps of our calculations 
we have ignored this kind of terms. This is not the case because this has been done 
in terms involving the b t functions, which contribute to the leading term, and mainly 
because the spin-dependent part is easily distinguished within experiments. 


In addition, for the hclicity amplitudes in this case we obtain 


^4++,++ 


0 + o(t ), 

4+o ,+0 


3 F u 

4+o,o+ 


j — r 9 F\ 03 

2 + 2 + - 

4+-,+_ 


0 + o(t ), 

4 +_,00 


/—7 i-9 F\ 03 

v t[ 2 2 

4+-,-+ 


0 , 

4o+,o+ 


463 

r 

~3x ’ 

4qo,oo 


_ /1 \ n 863 


3 (x ~ X ) Fl _ 34 


4.2 Range B: e ^«x« 1/VX 


yf-t (j/ ’i, 


yf—t 3 F\ , 


(71) 


The results we have obtained in our previous paper [3] in the small -x (but not exponentially 
small) regime for the polarized vector meson structure functions can be summarized in the 

14 The minus sign accounts for the difference in the four-dinrensional Minkowski metric components 7700 
and y zz . 
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following equations 


12 x 2 ^ 

F - 1 (T +71 
> *2 — -X~ fA + 1q) , 

bx 


L = 4^ h 

, h = —(A + J 0 ), 

(72) 

~{h + I o) 

c 

j A = — (A + A), 

bx 



9i ~ 0 > 92- ^(A + A), 


where A and A are ^-independent factors given by integrals of Bessel functions and some 
constants, plus the ^-dependence. According to these results, in this regime the Callan- 
Gross type relations become 


F -2 = 2 xF 1 (l + j-^j 


b 2 


2xbi 




(73) 


Notice that the first two relations have a Bjorken parameter factor x with respect to the 
ones in the regime A described in the previous section. The pre-factors are given by 


A 

A 


^ a 2 T, 2 q 4 (qR \ (2n+1)_D t 

VAA 4 T 2 (n + l)5 v’bJ n+1,D 

a 2 BY?q 2 /gi?\ (2n+3)_D 

TXA^f^Tl)W n ’ D ’ 


(74) 

(75) 


where C is a normalization constant and D — 2n + 1 + -^[2 A + a + (f + l)(p — 5)]. We 
have used the definition £ 2 = A 2 + A(A — 2 ). This is where the ^-dependence enters, since 
A = 7 (£)B — A. Note that the x-dependence of the structure functions is the same for all 
models, which was not the case in the previous section. The Bessel function integrals are 
defined as 15 


dujuj m K 2 n (uj) = /, 


n,m 


= 2 ? 


,_ 2 r> + n) i>- n) r» 1 

T( 2u) 2 v ’ 


(76) 


In the D4D8D8-brane and D4D6D6-brane models K n (w ) functions turn out to be fractional 
Bessel functions, however all the integrals are finite. Again, all these factors are independent 
of x and y 16 . Finally, the only model-dependent and ^-dependent numerical factor that 
enters the Callan-Gross type relations is 


5 = 1 + ^ = 1 + -^ 
A A+l ,D 


2D 


D + 2n + 1 


(77) 


15 Recall that this result is valid provided that the arguments of the gamma functions are positive. 
16 We refer the reader to [3] for a careful derivation and analysis of these results. 
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Model 

D3D7 

D4D8D8 

D4D6D6 

5 

2A+3 

A+2 

4(4A+3) 

8A+9 

16(A+1) 

8A+11 


Table 2: Parameter 5 in each model. 


Table 2 shows the values of 6 for the three models we consider. In all the cases we have 
1 < 5 < 2 since A is positive. As for the scalar mesons the results are the following 


F\ 

F 2 


7T 2 I |2^ A V 1 T 

16x 2p3 ' Ci ' U 2 / VF^N h,2A+3 ’ 

7T 2 2 /A 2 \ A_1 1 

~8x^ 3 Ci \?) °’ 2A+3+ i ’ 2A+s) - 


(78) 

(79) 


Equations (74) and (75) show the dependence of our small-:c structure functions with the 
momentum of the virtual photon: they always fall-off as a power of q 2 . In the Regge regime 
the calculations and measurements suggest that a Pomeron is exchanged 17 both at weak 
and strong ’t Hooft coupling. However, there is an important difference: when A is small 
the virtual photon scatters off a parton inside the hadron, leading to a small (growing) 
dependence in q 2 , while at large A a Pomeron is supposed to collide with the entire hadron. 
This is what happens in the present case, since we are using the gauge/string duality at 
strong coupling. Thus, we have obtained the same g 2 -power fall-off as in the large-:r case, 
and the amplitude is again dominated by the same double-trace operators [7]. 

Plugging our structure functions in the most general DIS differential cross section for a 
spin-one target we obtain 


da 


dxdydcj) 


MEe 4 
47 r 2 q 4 
MEe 4 


F-\ 

47r 2 g 4 
MEe 4 h 
487r 2 g 4 x 


(qq)[xy 2 + 2x(l - y)S\F 1 + 2xy 2 g 2 
(qq) [xy 2 + 2x(l — y)8} + 3 xy 2 5 
(qq)[y 2 + 2(1 - y)8] + 3y 2 5 


5- 

■si) 

i 

■ Sh){p 

•si) 

i 

■q) 


(p 

■q)‘ 

) 

(sh ■ 

ill _ 

(q- 

Sh)(p- 

si)' 

{.v■ 

q) 


ip- 

q) 2 


■si) 

_ (i. 

■ Sh)(p■ 

si)' 

i 

■q) 


(p- 

q) 2 


f 


(80) 


Once again, the experimentally relevant situation occurs when the lepton is longitudinally 
polarized, i.e. when si = Hik with Hi = ±1. Therefore, the differential cross section takes 
the form 

da MEe 4 Ii f. . r 2 n ,~ \ rl 35 rr , .1 , . 

{ te>[9 + 2(1 - ^ " me h ‘ s ^~‘ (81) 

17 There are some important differences between hard and soft Pomerons that we shall not address here, 
for this see [9]. 
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Then, we can study the behavior of this result for each hadronic polarization as follows. 


• Unpolarized target: ( qq) U n P oi = 1 , Sh = 0. In this case the spin-dependent term 
vanishes automatically and we obtain 


da ! 

"j 7 771 unpol 

dxdydcp 


MEe 4 /, 
487 r 2 g 4 x 


[: v 2 + 2(1 - i/M ' 


(82) 


• Longitudinally polarized target: (qq)LP = 0, = (0 ,H h Mz) 


da 

dxdydcj) 


LP 


0. 


(83) 


This is just like in the large-x regime. We obtain a vanishing result at leading order 
since ()\ — 0 (notice that the g x structure function together with the difference F x — \b\ 
give the main contribution to the longitudinally polarized cross section), and this is 
a consistency check of our results at small x for the polarized vector meson structure 
functions. 


• Transversally polarized target: ( qq)rp = 3/2, Sh = MH h x. Now, we have a vanishing 
Sh ■ k but a non-vanishing s-q — MHhk' x = ME'Hh sin 6 cos 0, and since E' = E(l — y) 
and sin(d) rsx Q rsj M/\q\ we finally obtain 


da 


MEe 4 /, f 3 


dxdydcj) 


\tp — 


48n 2 q 4 x 


M 


„[y + 2 (! - y)S\ - 3(5(1 - y)— HiHhCosct) 


(84) 


Since M/\q\ ~ \[M the second term also becomes sub-leading in this case. 


We can also obtain the helieity amplitudes in this regime which become 


t 4 ++i++ ~ 0 
71+0,0+ - 0 
7l+_,oo = 0{\/M) ^ 0 

7lo+,o+ = 2(5 — l)Fi > 0 


^4+o,+o — 3T 1 !, 

- 4 +-,+- - 0 , 

^4+-,-+ = 0) 

^00,00 = 3(5 — l)Fi > 0 


(85) 


Note that in this parametric region the factor <23 = 62/3 — 63 defined in the section 3.3 
becomes 03 = b 2 at leading order. This is interesting since it means that the b 2 structure 

function characterizes the double helieity flip amplitude. The corresponding _,_|_ helieity 

amplitude is proportional to t, thus it is always small in the DIS regime, but the b 2 ~ x~ 4 
growth could in principle lead to a non-vanishing helieity amplitude. 
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4.3 Range C: x ~ e ^ 

This is the exponentially small region of the Bjorken parameter. In the region B we have 
neglected the factor s 0 ' 1 ! 2 . The effect of this factor is very important within the parametric 
region x ~ exp (— y/X) [3]. The tildes indicate ten-dimensional Mandelstam variables. Let us 
consider the strong coupling regime 1 <C A <C N 18 , and exponentially large values of s with 
fixed. In this region the interaction cannot be considered local in the AdS background, 
since t becomes a differential operator [7, 9], and consequently we have to replace it by the 
following operator 

a't -> c/V 2 = a' , (86) 

which takes into account the momentum transfer in the transverse directions, and it induces 
the factor 

(a's) a ' i/2 ~ x~ a ' i/2 ~ x~ Q ' y2/2 . (87) 

The physical meaning of this is that within this region C there is a Pomeron exchange in 
the t-channel 19 . The consequence of this leads to include a factor x a '^l 2 multiplying the 
structure functions. This is obtained in the region C by replacing the differential operator 
—V 2 = — D 2 by its smallest eigenvalue (. Then, the structure functions become F\ oc 
x - 2 +a '(/2 an( ] Qt x~ 1+a '^ 2 , where the Pomeron intercept is identified as 2 — a'C t /2. This is 
the expected behavior for these functions. In particular, for the D3D7-brane model we have 
( = -^ 2(1 — A 2 ) < 0. This implies that in our results of the previous section 4.2 in the region 
B we have to multiply all the structure functions by the factor x a ’^ 2 . This implies that the 
differential cross sections for different meson polarizations as well as the helicity amplitudes 
must also be multiplied by this factor in the parametric region C. 

5 Results and discussion 

In this section we present our results and discuss them. 

5.1 Structure functions and differential cross sections 

Let us begin with the case of vector mesons. Figures 3, 4 and 5 show the structure functions 
F \, F- 2 , bi and b 2 for polarized vector mesons as a function of the Bjorken parameter. We 
consider the D3D7-, D4D8D8- and D4D6D6-brane models, respectively. In each case we 
focus on the lightest particles given by t = 1,2,3. For the small-x region the structure 
functions Fi, F 2 , b\ and b 2 are represented by a continuously decreasing curve, which has 

18 In this subsection we focus on the D3D7-brane model, the other models have a few modifications but 
the discussion is very similar. 

19 In reference [3] we have not considered multi-Pomeron exchange. 
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(c) bi(x) function (d) 62(0;) function 


Figure 3: Structure functions F\ , F 2 , b\ and 62 in terms of the Bjorken parameter x for the t = 1, 2 
and 3 vector mesons in the D3D7-brane model. Bell-shaped curves correspond to the large-x region 
(height decreases as l increases) while the hyperbolic curves correspond to the small-x region. 

very little dependence on the values of l = 1,2,3. On the other hand, for large x -values 
all the figures show a set of bell-shaped curves with maxima around x ~ 0.6. The height 
of these curves decreases as i increases. The normalization constants are specified in the 
next section and, they have been chosen in order to have the best fitting to the lattice QCD 
results of the first three moments of the p meson 20 . These moments are shown in table 5 and 
discussed below. Recall that in the small-x parameter region there is a model-dependent and 
^-dependent factor in the structure functions. All the structure functions behave similarly. 
On the one hand, the x > 0.2 region is dominated by a smooth bell-shaped function with 
maximum around x — 0.6 (the specific value depends on £, as it can be seen from equations 
(61) and (62)). It approaches zero as x —» 1. On the other hand, as x decreases these bell- 
shaped functions also decrease, and the small -x structure functions defined in equations (72), 
(74) and (75) dominate. The x < 0.1 region is described by functions x~ l or x~ 2 (depending 
on which structure function one considers). The two type of functions (the hyperbolic and 
the bell-shaped ones) overlap for 0.1 < x < 0.2, which is a transition region between both 
regimes of the Bjorken parameter. This behavior is understood from the Callan-Gross type 

20 We omit the possible dependence of these constants with l since we compare with the moments of the 
structure functions for £ = 1, i.e. the moments corresponding to the p meson, so the curves for t > 1 have 
been drawn for the same values of the constants as for £ = 1 . 
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(c) bi(x) function (d) 62(0;) function 


Figure 4: Structure functions F \, F 2 , b\ and b 2 in terms of the Bjorken parameter x for the t = 1, 2 
and 3 vector mesons in the D4D8D8-brane model. Bell-shaped curves correspond to the large-x 
region (height decreases as t increases) while the hyperbolic curves correspond to the small -x region. 

relations in each range of the Bjorken parameter, as well as from the fact that bi = 3F, 
in both cases. Also, remember that when x becomes exponentially small this description 
fails because the approximation that strings interact locally enough in order to use a flat 
space amplitudes breaks down. In this case the correct behavior should be obtained by 
considering a diffusion operator, and ultimately the description of full string theory on a 
curved background. The analysis in this parameter region is very difficult. However, there 
is evidence that the divergence disappears (at least for the glueball DIS [7]), leading to finite 
moments of the structure functions, but this is beyond the scope of this work. In the next 
section we will deal with this by introducing a low-x cutoff. We can see that for £ = 1 it 
gives similar curves for the three models we consider, even if the normalization constants are 
different for each model. It is interesting to see that the structure functions at small x have 
similar behavior for £ — 1,2, 3, while at large x they get smaller as the parameter £ increases. 
This is more evident in the D4D8D8 and D4D6D6-brane models than in the D3D7-brane 
model. In those cases, for i — 3 the corresponding curve is one order of magnitude smaller 
than the corresponding one for £ = 1 of the large-x structure function. 

We can also compare with the existing literature. In particular, at large x the curves for F\ 
and F 2 have a similar bell-shaped behavior as those presented in [12], although the maxima 
in that reference occurs for x ~ 0.8. In [12] the DIS structure functions for the lightest 


27 




















(c) bi(x) function (d) 62(®) function 


Figure 5: Structure functions F\ , F 2 , b\ and b 2 in terms of the Bjorken parameter x for the t = 1, 2 
and 3 vector mesons in the D4D6D6-brane model. Bell-shaped curves correspond to the large-x 
region (height decreases as t increases) while the hyperbolic curves correspond to the small- x region. 

vector mesons (they only consider the unpolarized case) have been obtained by using the 
D4D8D8-brane model, i.e. the p meson and the a\ axial-vector meson. In fact, in that 
reference F\ and F 2 structure functions have been obtained in the range 0.2 < x < 1. We 
recall that the calculations in that reference are different from ours since they considered a 
four-dimensional effective Lagrangian only describing the l/VX <C x < 1 range. Obviously, 
this describes only a part of the problem. Instead, in [1, 2] we have derived the corresponding 
Lagrangians for the three holographic dual models from first principles from type IIA and 
type IIB super gravities within this parametric region, while for small x in our paper [3] we 
have obtained the structure functions by using type IIA and type IIB string theories. From 
the structure functions obtained in [12] they found the ratio F 2 /(2Fi) « 1, satisfying this 
Callan-Gross type relation (as expected from the supergravity calculation [7] within this 
parametric range, note that it has not the factor x multiplying F\ which agrees with our 
results) within the interval 0.4 < x < 0.6 and momentum transfer 10 GeV 2 < \q 2 \ < 80 GeV 2 . 
However, in [1, 2, 3] we have proved that a modified version of the Callan-Gross relation 
holds within the whole range of the Bjorken parameter and for any value of the momentum 
transfer, in fact there appear two Callan-Gross type relations: F 2 = 2F\ for l/\/A <C x < 1 
(supergravity) and F 2 = 2xF 1 6 for exp(—\/A) <C x <C l/\/A (from superstring theory, see 
section 4.2). Note that 0.5 < 5 < 1 for all mesons described by the models we consider. 
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Moreover, we have also obtained new exact relations among the structure functions for scalar 
and for polarized vector mesons in the planar limit at strong’t Hooft coupling. In addition, 
in reference [12] only structure functions for the lowest states of vector and axial vector 
mesons were studied, while our results hold for any scalar or polarized vector meson. 

Now, we present our results for the DIS differential cross section of unpolarized vector 
mesons depicted as contour line maps in terms of x and y variables. We show that for large 
x in figure 6 and for small x in figure 7, and for the three different Dp-brane models we have 
considered. Broader curves in figure 6 correspond to lower values of the DIS differential 
cross section. On the other hand, in figure 7 curves with larger slope correspond to larger 
values of the DIS differential cross section. In figure 6 we can observe that the differential 
cross section of unpolarized holographic vector mesons is larger around x ~ 0.5, and also it 
becomes less spread for the Sakai-Sugimoto model. Another feature is that this differential 
cross section is larger for smaller values of the fractional energy loss of the charged lepton 
y , having a maximum when the energy of the incident lepton and the one of the scattered 
lepton are the same. On the other hand, curves in figure 7 show contour line maps of DIS 
differential cross section for unpolarized holographic vector mesons for the three flavor-brane 
models we describe, in the region exp (—\/A) l/VX The horizontal axis represents 

the Bjorken parameter within the range [0.01, 0.1], while the vertical one corresponds to the 
variable y, [0, 1]. 

Next, let us focus on the comparison with phenomenology of the meson structure functions. 
This turns out to be difficult because the experimental data of these functions is actually 
much more limited than in the case of baryons. Moreover, it is more difficult to carry out 
such a comparison when the Bjorken parameter x is small. In the case of pions there have 
been important phenomenological studies [13, 14, 15, 16], however all these investigations 
focus on the valence structure function related to fixed-target pion Drell-Yan experiments. 
The corresponding data are restricted to the parametric region x > 0.2. Notice that in 
our calculations the pion corresponds to the l = 1 meson in the D3D7-model, since it is 
the lightest pseudoscalar particle. These phenomenological studies could in principle be 
compared with the results obtained in our previous papers studying the parametric range 
lfy/X <C x < 1 [1, 2], The small- x behavior we find is difficult to compare with the existing 
phenomenological models. However, in [17] it was proposed a model in which the parton 
distribution functions have a sea term that diverges as x —> 0, as it occurs for our structure 
functions. In addition, this type of divergences beyond the valence term was shown to appear 
on the pion contribution to the deuteron structure function b\ (see for instance figure 3 in 
reference [18]). 

In general terms, our approach shares fundamental concepts with references [9, 19, 20, 
21, 22] as well as it involves somehow similar techniques to the ones used by them. Very 
interestingly, in those references there is agreement with baryon DIS data from HI and 
ZEUS experiments at HERA. In fact in these references the objects of study are the glueball 
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(a) D3D7 & £ = 1 (b) D4D6D6 & £ = 1 (c) D4D8D8 & £ = 1 



(d) D3D7 & £ = 2 (e) D4D6D6 & l = 2 (f) D4D8D8 & £ = 2 



(g) D3D7 & £ = 3 (h) D4D6D6 & £ = 3 (i) D4D8D8 & £ = 3 



(j) D3D7 & £ = 10 (k) D4D6D6 & £ = 10 (1) D4D8D8 & £ = 10 


Figure 6: Contour line maps of DIS differential cross section for unpolarized holographic vector 
mesons for the three flavor-brane models we consider, in the region l/\/A <C x < 1. Different 
mesons are labeled with £ = 1,2,3 and 10 (displayed on the four rows). The differential cross 
sections have been normalized as explained in the main text. The horizontal axis represents the 
Bjorken parameter within the range [0.1, 1], while the vertical one corresponds to the variable y, 
[0, 1]. Broader curves correspond to lower values of the DIS cross section. 
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(j) D3D7 & t = 10 


(k) D4D6D6 & £ = 10 (1) D4D8D8 & l = 10 


Figure 7: Contour line maps of DIS differential cross section for unpolarized holographic vector 
mesons for the three flavor-brane models we consider, in the region exp (—\/A) <C x <C l/VX 
Different mesons are labeled with i = 1,2,3 and 10 (displayed on the four rows). The differential 
cross sections have been normalized as explained in the main text. The horizontal axis represents 
the Bjorken parameter within the range [0.01, 0.1], while the vertical one corresponds to the variable 
y , [0, 1]. Curves with larger slope correspond to larger values of the DIS cross section. 
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structure functions. Since the structure functions of scalar mesons behave similarly as those 
of glueballs for all values of the Bjorken parameter, we regard it as an indication suggesting 
that our predictions should be correct. Keeping in mind that the results for glueballs can be 
compared with experimental information from HERA, we would expect that our results for 
mesons (which at least for scalar ones are closely related to the glueball results of [7]) could 
be confirmed. 


5.2 Comparison with lattice QCD 

Now, let us focus on the comparison with lattice QCD results. Particularly, we can compare 
the lower moments of D1S structure functions of the pion and the p meson obtained from 
lattice QCD [23] with our own predictions. The n-th moment of a generic function h(x, q 2 ) 
is given by 


M n [h]= / dxx n 1 h(x,q) 


( 88 ) 


We consider the spin-dependent structure functions of the pion and p mesons. In particular, 
the lattice QCD calculations in [23] have considered Wilson fermions and they have been 
carried out for three values of the quark mass. Therefore, it allows one to carry out an 
extrapolation to the chiral limit. The way we can compare our results with the lattice QCD 
ones for the pion is by using our expressions of the F 2 structure function 21 for small and 
large x values, and integrate them with appropriate lower and upper cut-offs. The former is 
integrated from x = 0.0001 to x — 0.1, while the latter is integrated from x = 0.1 to x — 1. 
Notice that within each structure function there is an undetermined constant containing A, 
R , q 2 and some other numbers. In fact, there is a different constant for each parametric 
regime of the Bjorken variable. Thus, for each model there are two undetermined constants. 
They are effectively two free parameters for each flavor-brane model, which are determined 
by carrying out the best fitting to the lattice QCD data for the moments Mj(F 2 ) with 
i = 1,2, 3. The rest of the fittings are performed in the same way. 

Our results and the lattice QCD data as well as the deviations from each other are listed in 
table 3. The constants of the small -x regime which leads to the best fitting to the lattice QCD 
data of reference [23] are 0.025, 0.020 and 0.022, while for the large-x regime we found 24, 
2020 and 267, for the D3D7-, the D4D8D8- and the D4D6D6-brane models, respectively. We 
find discrepancies with respect to the lattice QCD calculations up to 18% for our results in 
the case of the pion, £ = 1. Recall that in the 1/N expansion discrepancies are of order 30%, 
thus our results are within the expectations. We should keep in mind that when analyzing 
discrepancies between gauge/string duality predictions, phenomenology and experimental 
data, for instance bottom-up models in Eve dimensions (the AdS/QCD model) give results 


21 For the pion we use F 2 and not F\ because one of the features of the large N limit is that the virtual 
photon strikes the entire hadron, thus F\ = 0 in the large-a; region. 
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Model / Moment 

M!(F 2 ) 

M 2 (F 2 ) 

M 3 (F 2 ) 

Lattice QCD 

0.3047 

0.1180 

0.0583 

D3D7 

0.3067 

0.0962 

0.0658 

Percentage error 

-0.6 

18.5 

-12.8 

D4D8D8 

0.3061 

0.1018 

0.0643 

Percentage error 

-0.5 

13.7 

-10.3 

D4D6D6 

0.3064 

0.0990 

0.0650 

Percentage error 

-0.6 

16.1 

-11.6 


Table 3: Comparison of our results for the first moments of the structure function F 2 of the pion 
for a suitable choice of the normalization constants with respect to the average results of the lattice 
QCD computations in [23]. Uncertainties in the lattice computations are omitted. 


related to masses and decay constants of mesons with discrepancies of order 5% [24, 25]. 
On the other hand, for the case of more involved calculations leading to the A/ = 1/2 rule 
describing the kaon decays, related to the calculation of four-point correlation functions, 
discrepancies become of order 30% [26, 27]. 

There are more recent results in the lattice QCD literature to compare with. For instance 
we can consider data from references [28, 29]. The results for these fittings are shown in 
table 4. One can see that all the results are improved being the deviation smaller than 10%. 
In this case, the constants we find are 0.014, 0.011 and 0.012 for the small-re regime, and 29, 
2450 and 325 for the large-re regime, for the D3D7-, the D4D8D8- and the D4D6D6-brane 
models, respectively. 

It is possible to understand a bit better the comparison between lattice QCD data and 
our holographic results by isolating the contributions from the small -x and large -x struc¬ 
ture functions to each moment. We are tempted to associate the small-x region with the 
concept of sea-quark distribution function, which describes the possibility of finding a quark 
or an anti-quark (generated by gluon splitting) carrying a very small fraction of the hadron 
momentum. In this picture the large -x structure function is then related to the valence 
distribution function associated with quarks that carry a considerable fraction of the hadron 
momentum 22 . Table I in [30] shows that both for lattice and for phenomenological results 
one finds that the contribution of the sea of quarks is important only for the first moment 
of the structure functions, since it yields a considerable fraction of the final result. In the 

22 These are not formal definitions of these concepts. 
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Model / Moment 

Mi{F 2 ) 

M 2 (F 2 ) 

M 3 (F 2 ) 

Lattice QCD 

0.27 

0.13 

0.074 

D3D7 

0.2708 

0.1161 

0.0803 

Percentage error 

-0.3 

10.7 

-8.5 

D4D8D8 

0.2705 

0.1221 

0.0779 

Percentage error 

-0.2 

6.1 

-5.2 

D4D6D6 

0.2706 

0.1191 

0.0791 

Percentage error 

-0.2 

8.4 

-6.9 


Table 4: Comparison of our results for the first moments of the structure function F 2 for the pion 
for a suitable choice of the normalization constants with respect to the average results of the lattice 
QCD computations in [28, 29]. Uncertainties in the lattice computations are omitted. 


second and third moment its contribution is substantially reduced, and the contribution from 
the valence distribution function gives almost the full result. This is similar to what happens 
with our structure functions: the integral for the small-x region is important in order to fit 
the first moment Mi, but one can almost ignore it for M 2 and M 3 . 

Table 5 shows the comparison of our results for the first moments of the structure function 
Fi for the p meson for a suitable choice of the normalization constants with respect to the 
average results of the lattice QCD computations in [23]. The constants in this case are 0.013, 
0.011 and 0.012 for the small-x regime, and 14, 1229 and 162 for the large-x regime, for the 
D3D7-, the D4D8D8- and the D4D6D6-brane models, respectively. As it was commented 
for table 3, it is very interesting that the relative discrepancies are smaller than 21 percent. 
Now, by using tables 4 and 5 we can draw figures 8 and 9. In figure 8 the first three moments 
for F 2 are displayed for the pion. The free parameters for each dual holographic model are 
chosen in order to fit the results in [23] obtained with lattice QCD as explained before. The 
three studied holographic dual models are shown, together with the value obtained in the 
reference. We find good agreement with lattice QCD data with discrepancies lower than 
20%. By using more recent lattice QCD results of references [28, 29] the fittings improve 
with discrepancies lower than 10%. In addition, in figure 9 the first three moments of F\ 
are shown for the p meson. The three holographic dual models are presented in a similar 
way as in figure 8. Discrepancies with respect to that reference are under 21% in all cases. 
Therefore, this also shows good agreement with lattice QCD data. 

Notice that in all our fittings the constant of the structure function corresponding to the 
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Model / Moment 

M 2 (F 1 ) 

M 3 (F i) 

M 4 (Fi) 

Lattice QCD 

0.1743 

0.074 

0.035 

D3D7 

0.1755 

0.059 

0.040 

Percentage error 

-0.7 

21.3 

-14.1 

D4D8D8 

0.1752 

0.062 

0.040 

Percentage error 

-0.5 

16.4 

-11.8 

D4D6D6 

0.1754 

0.060 

0.040 

Percentage error 

-0.6 

18.8 

-13.0 


Table 5: Comparison of our results for the first moments of the structure function F\ for the p 
meson for a suitable choice of the normalization constants with respect to the average results of 
the lattice QCD computations in [23]. Uncertainties in the lattice computations are omitted. 


F 2 n moments 


7r-meson moments 


0.25 


0.20 


0.15 


0.10 - 

0.5 1.0 


• D3D7 
■ D4D8 

♦ D4D6 
▲ Lattice 


f 

1.5 2.0 2.5 3.0 


— 1 n-th moment 
3.5 


Figure 8: The first three moments of F 2 are shown for the n meson. The free parameters for each 
dual holographic model are chosen in order to fit the results in [23], obtained with lattice QCD 
as explained before. The three holographic dual models are shown, along with the average value 
obtained in the mentioned reference. Discrepancies with respect to the reference are under 20%. 
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p-meson moments 


F\ p moments 

0.18 r 
0.16 1 
0.14 1 
0.12 1 
0.10 1 
0.08 1 
0.06 1 
0.04 1 

0.5 1.0 1.5 


C 

. ? 

2.0 2.5 3.0 


• D3D7 
■ D4D8 

♦ D4D6 
▲ Lattice 

— 1 n-th moment 
3.5 


Figure 9: The first three moments of F\ are shown for the p meson. The free parameters for each 
holographic dual model are chosen in order to fit the results in [23], obtained with lattice QCD. The 
three holographic dual models are shown, along with the average value obtained in this reference. 
Discrepancies with respect to the reference are under 21% in all cases. 


small -x regime is always a few orders of magnitude smaller than the one associated with the 
large-a: structure function. This is in full agreement with our calculations, at least in the 
regime where the holographic dual model represents an accurate description of the dual gauge 
theory. This could have been directly predicted from our analytical results. By examining 
the difference between the constants in front of (65) and (78, 79) one can see that for the 
£ = 1 scalar meson the numerical coefficients and normalization constants are of the same 
order, and the main difference is given by the A~ 1//2 factor that appears in the small -x case. 
This is small at large ’t ffooft coupling, and it is directly related to the value of a: = 0.1 
which we use in order to separate the first two regimes (A and B) in our calculations. 

Another important point comes also from the comparison between our analytical results 
and the different models formulated in the literature in order to fit the phenomenological 
data for the measurement of the pion valence distribution functions. As it is explained 
in [31, 30] there have been different attempts to describe the available data, and while 
there are some differences between them, they all agree on the fact that the distribution 
functions should behave as (1 — x) 2 in the x —> 1 limit. In [32] it has been found a fall-off 
(~ (1 — x) 2±01 ) for the valence distribution function consistent with Drell-Yan data. This 
is consistent with theoretical predictions based on perturbative QCD (see for instance [33] 
and references therein) as well as with calculations using Dyson-Schwinger equations [34], 
This is exactly what happens with our structure functions (65) in the D3D7-brane model for 
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the i = 1 scalar meson, i.e. the pion. This would suggest that the D3D7-brane model gives 
better predictions than the other two ones studied in this work. In the D4D8D8-brane and 
D4D6D6-brane models the exponents are 4.59 and 3.33, respectively. 

A very interesting review of lattice QCD calculations on moments of hadron distribution 
functions, as well as other important hadron structure observables is given in [35]. The 
earliest calculations of the lowest two moments of the pion quark distribution function have 
been done in [36, 37, 38] in quenched lattice QCD. An extensive study of several moments of 
the pion and the p meson have been done in [23]. In fact, we have compared our results with 
those in that reference in our tables 3 and 5 and in our figures 8 and 9, while in our table 
4 we compare with references [28, 29]. The earliest results of [23] have been obtained by 
using quenched lattice QCD simulations, which means that the fermion determinant is set 
to one, thus neglecting all quark-antiquark quantum fluctuations. This introduces an error 
which can only be quantified by comparison with dynamical simulations with physical quark 
masses, which are only computationally possible since a few years 23 . In fact for instance in 
[28, 31, 29] dynamical quarks have been considered. 

We may also try to use this kind of fittings to compare the first few moments corresponding 
to our bi functions with the ones calculated in [23]. However, in that paper the results found 
for these moments are smaller than the Fi moments 24 . Since the structure functions we 
obtained from string theory satisfy the relations bj = 3iq, the numerical comparison becomes 
meaningless. It seems that there is a conceptual difference, and it would be very interesting 
to study where it comes from or to be able to resolve this discrepancy from new experimental 
data. 

Once again we would like to emphasize that the discussion presented in this last section 
has to be understood qualitatively. One should keep in mind that the gauge theories that we 
are able to study in terms of string theory dual models do not look exactly like QCD. What 
makes them interesting though is that they have some important similarities related to the 
meson spectrum and the concepts of color confinement, chiral symmetry breaking and other 
properties. Thus, it is really interesting that these holographic dual models give predictions 
of the moments of the structure functions which are comparable to lattice QCD results. 

The good agreement we have obtained with lattice QCD data is somehow not unexpected 
taking into account several examples where the supergravity dual models give results which 
compare reasonably well with lattice QCD data. For instance, in the case of SYM plasmas 
good level of agreement has been found between lattice QCD results and the gauge/gravity 
duality approach by including string theory corrections to the supergravity calculations as 
it has been shown for the mass and electric charge transport coefficients. In the case of the 
mass transport the quantitative agreement between lattice QFT calculations for rj/s [39] and 

23 We thank Andreas Schafer for this comment. 

24 In fact, even if we actually tried to do this the very small values of these moments and the considerable 
errors associated with them in [23] would probably make the fittings unreliable. 
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the gauge/string duality [40, 41] is remarkable (r/ is the J\f = 4 SYM plasma shear viscosity 
and s the entropy density). Also for electric charge transport the DC conductivity has 
been calculated using lattice QCD [42] and the gauge/string duality [43], including 0(a' 3 ) 
corrections, leading to good level of agreement [44, 45, 46, 47]. Of course, also in those 
cases the gauge/string dual calculations strictly hold for large N and the strong coupling 
expansion (1/A), but if one insists on using N — 3 and A ~ 15 it leads to the results we 
mention in this paragraph. 

Also, the fact that the methods developed for the holographic dual description of DIS can 
be extended to other strongly coupled SU(N) gauge theories is very interesting. 
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